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We demonstrate that for the sine-Gordon theory at the free fermion point, the 2- 
point correlation functions of the fields exp(ict$) for < a < 1 can be parameterized in 
terms of a solution to a sinh-Gordon-like equation. This result is derived by summing 
over intermediate multiparticle states and using the form factors to express this as a 
Fredholm determinant. The proof of the differential equations relies on a Z2 graded 
multiplication law satisfied by the integral operators of the Fredholm determinant. Using 
this methodology, we give a new proof of the differential equations which govern the spin 
and disorder field correlators in the Ising model. 
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1. Introduction 

It is well-known that for certain integrable quantum field theories in 2 dimensions, 
non-trivial correlation functions can be parameterized by solutions of non-linear partial 
differential equations^. The first result of this nature was obtained by Wu et. al.0, where 
it was shown that the 2 point correlation functions of spin and disorder fields can be 
parameterized in terms of a solution to the Painleve III equation. Subsequently, this result 
was obtained in a variety of ways @@]. Miwa, Jimbo and Sato obtained the result starting 
from the braiding relations of the spin fields with the free fermion fields. In 0], using the 
form factors, the sum over intermediate multiparticle states was related to tau-functions 
satisfying known differential equations. The two latter derivations are not unrelated, since 
in || the form-factors are implicit in the expressions for the spin fields in terms of the 
fermions. 

Differential equations for correlation functions in other models which possess free 
fermion points were derived by Its et. al.00. In this approach, the Bethe ansatz is used 
to express the correlation function as a Fredholm determinant, and general techniques were 
developed for obtaining differential equations from this representation. 

In this paper, we derive differential equations that characterize the non-trivial corre- 
lation functions of the sine- Gordon theory at the free fermion point (SGff). The method 
we use is a kind of hybrid of the methods discussed above. Namely, using the form fac- 
tors of the fields, we express the sum over intermediate multiparticle states as a Fredholm 
determinant. The type of kernels we obtain differ from the ones considered in [H[BI. Nev- 
ertheless, we find that the methods developed there can be generalized to deal with the 
relativistic kernels which generically arise from the form factor sum. Indeed, using this 
methodology, we also give a new proof of the differential equations for the Ising theory. 
In should be pointed out that Fredholm determinant expressions in the SG / f theory were 
also obtained in , however the latter determinants have a different origin than the form 
factor sum, namely they originate from the Bethe ansatz, and do not coincide with our 
expressions. 

3 This manuscript is a corrected version of our work which has appeared in Nucl. Phys. B426 
(1994) 534. Previously, we had not found the most general solution to the differential equations 
(3.37), but rather projected onto a solution which is valid only for a = a' . The only changes in 
this new version are at the end of section 3, where the general solution for a ^ a' are obtained. 
We thank S. Lukyanov for drawing our attention to a possible error. Recently, similar results 
were obtained by H. Widom using different methods]!]]. 
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In || the correlation functions of fields with non-trivial monodromy with charged 
free fermion fields (as in eq. 2.14 below) were studied. Since such fields are essentially 
completely characterized by this monodromy, one can interpret the work || as dealing 
with the sine-Gordon theory, and the results on differential equations presented below are 
implicit in ||. The methods we use are different and perhaps more direct. Generalization 
of the methods employed in [Q to the n-point functions is discussed in 

The SGff theory is perhaps the simplest relativistic quantum field theory where all 
of the integrability structures have been developed, such as the Lax pair, classical and 
quantum monodromy matrix, quantum inverse scattering, etc., and furthermore has non- 
trivial correlation functions of the basic fields on which these structures are based. This 
is to be contrasted with the situation in Ising, where even a lagrangian involving the spin 
fields is unavailable. For this reason, we believe our SGff results may be instrumental for 
investigations beyond the free fermion point. 

The result we derive is simply described. The sine-Gordon theory is described by the 
Euclidean action 

S = J d 2 z (d z $dj$ + A cos(^$)) , (1.1) 



where z = (t + ix)/2, z = (t — ix)/2. The free fermion point occurs at /3 = 1 JT^] . Consider 
the correlation functions 

(0| : e ia *(*>*) : : e *«'*(0) : | ) = exp E(z, z), (1.2) 

where < a, a' < 1. Since the bosonization relations are : exp(±z$) := V'iV'ip where 
ip±, are the free massive Dirac fermions, the fields : exp(io;$) : for < a < 1 are not 
simply expressed in terms of the fermions, and thus have non-trivial correlation functions. 
We show that these correlation functions can be parameterized by a solution <p(r) of a 
differential equation 

d 2 r + -dr] E(r) = — (1 - cosh2y?) 
r I 2 



(1.3) 

_o J. _ \ lib . <±ICt — LX I , i 2 \ 

o* H — o r I cp = — - smh 2cp H ^ tanh99(l — tanh 99), 



1 „ \ m 2 . , „ 4(a - a' yJ 



r 



2 ^ r 2 

where m is the mass of the free particles, and r 2 = 4zz. We will also show that the 
result (|l~3|) continues to be true if one replaces the vacuua in (|1.2|) with multiparticle states. 
The equivalence of the differential equations ( |1.3|) and the Painleve ones is well known. 
(See e.g. [0.) 

As discussed below, when a, a' = ±1/2, the SGff correlation functions can be related 



to squares of Ising correlation functions [ 12 ] [Oj. Our results are obtained entirely within 



the SG f f theory and do not use this connection with Ising 2 . 
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2. Correlation Functions as Fredholm Determinants 



2.1 Form Factors 

The particle spectrum of the SGff theory consists of free solitons with topological 
U(l) charge e = ±1, where the solitons are identified as the free charged fermions. The 
N-particle states will be denoted as 

• • • , 0iv)ei- ej v> ( 2 -l) 

where 9{ (e*) is the rapidity (U(l) charge) of the i-th particle, E{ = mcosh6>i, Pi = 
msinh^. Inner products are normalized as follows: 

S '(9'\9) £ = 5 S £ ' 5(9-9'). (2.2) 

One then has the following resolution of the identity: 

00 -1 n OO n OO 

1= J2mJ2 / Mi'" d9 N \9 1 ---9 N ) ei ... eN e "- s *{9 N ---9 1 \. (2.3) 

Consider now the quantum fields cr a , 

a a (z,z) = : e iaHz ^ : (2.4) 

where $ is the SG field with the normalization implicit in the action ( |1 . 1| ) . These fields 
are local, Lorentz spinless, and U(l) neutral. They have the hermiticity property: 

(a a ) f =a_ Q . (2.5) 

In terms of the U(l) charged fermion fields ijj±,if)±, the bosonization relations are 
: exp(z$) := iJj+iJj_, : exp(— i$) := ipj^ij)-. Thus, we limit ourselves to the range < a < 1. 
Fields outside this range can be obtained through operator products with fermion bilinears: 

o~ a ±i ~: ip±ip T cr a ■ (2.6) 

The form factors of these fields are defined as 

F a (9 u ■ • • , 9 N ) £l ... £N = (O|<7 a (O)|0i ■ • ■ 9 N ) £1 ... £N . (2.7) 
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Since the fields a a are U(l) neutral, the form factors are non-vanishing only for U(l) 
neutral states. Explicit expressions for the form factors are as follows. For convenience 
define 

fa(@l, " " " 5 #2n) = F a (8i, • - • , 02n)e!— e 2 „> e i = — £ i+n = ~ lj i = 1, ••) (2.8) 

It is also helpful to define the momentum space variable 

u = exp(9). (2.9) 

Then one has 

/«(*!• "flan) = (Cam)" 2 (_l)»(»-D/a (^^f ^II(^)^ a (^ + n)^ a 



(2.10) 



x A(iti, • • ■ , Any- 
where 



a / v FL^XnK - U j) Un+l<i<j( U i ~ U i) , 01 - 

A(wi, • • ■ , An) = - 2n ~ (2.11) 

rir=llL=n+lK + A)- 

2 

The factor (c Q m) a is included based on our knowledge of the scaling dimension of 
the field <r a , which is the same as in the massless limit. The dimensionless constants c a 
can be fixed from knowledge of the 1-point functions: 

(0K(0)|0> = (c a m)- 2 . (2.12) 

We henceforth set c a m = 1. 

The form factors for other assignments of the charges e$ are trivially related to f a by 
the relation 

F a (9i, ■ ■ ■ , 8i, di+i, • • • , ^jv)ei---e i e i+ i---e iv = —F a {0\, • ' ' , 6i+i,9i, ■ ■ ■ , 0jv) £l ...e i+ ie i — ejv • 

(2.13) 

The above form factors can be derived in a variety of ways. They were first obtained 
implicitly by Schroer and Truong[|l3| by first expressing the SG field as a non-local fermion 



bilinear using d z <& = 9j$ = and then normal ordering the expressions 

exp(za$) with respect to the fermions. For a = ±1/2, these form factors were also ob- 
tained by Smirnov using the bootstrap axioms [0. In appendix A, we give an independent 



derivation of the form factors based on the monodromy relations: 

a a (x,t) ip±{y,t) = ip±(y,t) a a (x,t) for y < x 

(2.14) 

= e ±2ma ip±(y,t) a a (x,t) for y > x 
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The latter derivation is an extension of techniques developed by Jimbo, Miwa, and Sato 
for the Ising theory 0]. Finally, we mention that the form factors can also be derived using 



vertex operator techniques in radial quantization |15[ . 

We now describe some of the main properties of the above form factors. The her- 
miticity relation (|2.5f) necessarily implies 

F a {9\ • ■ ■ 6>Ar) ei ... ejv = F- a (0N ■ ■ ■ #i) e ^... e *, (2-15) 

where e* = —e. The expression ( |2.10| ) is primarily characterized by the residue property: 

2m f a {6 x ■ ■ ■ d n , n +z7T+ 7 , 9 n+2 ■ ■ ■ 9 2n ) ~ ^ ~ 6 ^ fa{0i • ■ ■ 0n+2, ■ ■ ■ (2-16) 

7 

The above relation is a generalization of the residue bootstrap axiom for form factors. 



As explained in |16|] , the additional phase exp(27ua) is a consequence of the monodromy 



2.2 Fredholm Determinants 

Consider first the 2-point correlation functions. Inserting the resolution of the identity 
one obtains 



oo ~ . / N 

(0\a a (z,z)a al (0)\0)= ]T — ^ / d9 1 ■ ■ ■ dd N J] exp(- 

N=0 ' £!,-••, £jv \i=l 



-mz Ui — mz u- " 



x F a (6i • • • 0N)ei—e N Fa'(9N " " " &l)e* N — e* 



? w / ^ " ' d02n ( n exp( 



-mz it,- — mz «„■ 1 



(2.17) 

In the second equation, we have set N = 2n, and summed over the P oss ^ e P er_ 

mutations of U(l) charges. 

We now show how these correlation functions may be expressed as a Fredholm deter- 
minant: 

(0\a a (z,z) t7 a /(0)|0> = Det(l + K (a ' a,) ). (2.18) 
Above, K^ a ' a ) is a 2 x 2 matrix of Fredholm integral operators: 

(£±3 £;tD- (2 - i9) 



For simplicity of notation, we will not always display the a, a' (and z, z ) dependence of the 
kernel K. The basic Fredholm theory of scalar integral operators can be easily extended 
to finite matrices of integral operators. One obtains 

' dui duN 



Det 1 + K 



j\f] ^2 

N=0 £i---£jv 



«1 



UN 



det{K eijej (u h Uj)}, 



(2.20) 



where det{K £u£ . (ui, %)} is an ordinary determinant of the finite N x N matrix with 



it,-, uA as entries: 



det {K ei>E:j (ui,Uj)} 



K eiei («i,iti) K £l£2 (u ll u 2 ) 

K £2£l (u 2 ,U 1 ) 



K 



El ,£JV 



[Ui, U N ) 
K £2 e N (u2,U N ) 

K Sn£n (un,un) 



(2.21) 



K £N£l {u N ,ui) 

With the appropriate choice of kernel K, the expression ( ^.2U ) coincides with the form 
factor sum (|2.17|) . The result is 

K ++ (u : v) = K—(u,v) = 
e(u)e(v) 



K- + {u,v) 
K+-(u,v) 



U + V 

e(u)e(v) 

U + V 



(2.22) 



where 



e(u) 



e(u) 



sin(7ra) 



7T 



1/2 



u (l+ a '-a)/2 exp ^_™ ( ^ + - u -l^ 



sin(7ro; / ) 



TV 



1/2 



(2.23) 



ti 



(l+a-a')/2 



/ m . _ 

exp I (zu + zu 

\ 2 



The above choice of kernel is not the unique one that reproduces the expression ( |2.17|) . 
However, viewing ( ^.22| ) as an ansatz, then the choice of e, e is unique up to constants. 
The result (|2.22| )( pT2"3|) is proven as follows. Note first that the condition K ++ = 

K = automatically projects the sum over £7(1) charges e% in ( |2.20|) onto neutral 

states. Let 



{K Si£ (ui,Uj)} 



so that 



A 
B 



Si Ei^fi 1, % 1, .., 77., 



(2.24) 



A 
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B 



K + _(v 2 , ui) 



\K + -(v n ,ui) 



K + -(v 1 ,U n ) 

K + -(v 2 , u n ) 



K + -(v n ,u n )J 



where we have defined Vi = Ui+ n . Using the Cauchy identity, 

1 



det 



A(lti ■■■U2r 



Ui + Vj 

we write the determinant of {K £i£j (m, Uj)} as 

det{K eiej (ui,Uj)} = det ( ° f) = (-l) n det A det B 

n 

= (-l) n A( Ul ■ ■ ■ u 2n ) 2 \\ (eiujeivi))' 



[2.25) 



(2.26) 



Comparing ( p.!7|) and (|2.20|) , and using the explicit expressions (|2.10|) of the form factors, 
one sees that to prove eq. ( |2.18| ) we need to satisfy the relation 



:-l)-J](e( Mi )e(^)) 2 = J] ((^">^ +c >^ +Q '(^) W 

x (exp(— mziui + u») - vriz{u~ x + w" 1 ))) 



i=l 



i=l 



(2.27) 



The functions e(u),e(u) given in ( |2.23|) were easily chosen to satisfy ( |2.27j ). 
Consider now a correlation function of the kind 

($V a (z,z)<7 <(0,0)|$), 



(2.28) 



where and (^'\ are two C/(l)-neutral eigenstates of the Hamiltonian; i.e they are pure 
multi-particle states = \^i,---,^N)e 1 ,---,s N with J2i £ i = 0' an d similarly for \^'). 
The correlation function ( |2.28|) can be computed using the crossing properties and the 
factorized expression ( |2.1U| ) of the form factors. One has : 



oo ^ „ / In 

(y'\a a (z, z)a a , (0,0)\V) = ^ — — / d9 l ---d9 2n \ l[exp 

n=0^ H -> J \i=l 



—mz Ui — mz tt„- 



2n 



(2.29) 



X 



fa(0i ■ ■■d 2n )f a >(d 2n •$!) J] fa^fr; 



where f a -,a'(Q', is a complicated function depending on the rapidities (3 and (3' . Since 

the extra dependence on \& and is factorized into the product Yli fa;a'(^i'i it can 

be absorbed into a redefinition of the functions e(u) and e(u). Thus one still has 

(*V«(*,;?K'(0,0)|W) = Det (1 + K a , a / ; *,*0 , (2-30) 



where K is still of the form ( |2.22 ), but with more complicated functions e, e. However, 



their (2, z) dependence are identical. This is the origin of the fact that all these correlations 
satisfy the same differential equation. 

2.3 Relation to Ising 2 

Since the Ising model is equivalent to a free real Majorana fermion, whereas the 
SGff theory is a free complex Dirac fermion, at the lagrangian level, the SGff theory is 
two copies of the Ising theory. It is well known that the above observation is somewhat 
naive for a number of reasons. For example, in the massless limit, the Dirac and Ising 2 
theories differ in their detailed partition functions on the torus [ff7|. For our purposes, the 



limitation of the identification of SGff with Ising 2 is that only correlation functions of 
certain special fields in the SGff theory are related to Ising 2 correlation functions. Let 



a(z, z) and \i{z,z) denote the Ising spin and disorder fields respectively. In [|12|JT3[ the 
following identifications were made: 

(0\cr(z, z)o (0)|0> 2 = (0|sin($0,:?)/2) sin($(0)/2)|0) 

2 ( 2 - 31 ) 

(0|M^)M0)|0> 2 = (0|cos($(z,Z)/2) cos($(0)/2)|0), 

where on the RHS are SGff correlation functions. In this section, we present a different 
derivation of part of the relations ( |2.31| ) based on Fredholm determinants. 



We first review some aspects of the form factor sum for the Ising correlators |18| [|19 
Define 

t±(z,z) = (0|M^)M0)|0> ± (0\a(z,z)a(0)\0). (2.32) 

In the free fermion description, the spectrum consists of a single neutral particle. From 
the known form factors of the spin and disorder fields, one obtains 



n=0 N ' i<j N J ' \i=l 



mz Ua 1 



(2.33) 
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The functions t± may be expressed as a Fredholm determinant for scalar integral operators: 



r± = Det(l±V r ), (2.34) 

where V is the kernel 



V(u,v) = ^-(^) expl-^izu + zu- 1 )). (2.35) 
27T \u + v I V 2 



In terms of the functions t±, the relations ( 2.31 ) read 

(0\a 1/2 (z,z)a 1/2 (0)\0) 



1 (2-36) 

(0|ffi /a (z,z)(7_ 1/a (0)|0> = -(7j + 7^) ) 

where we have used the hermiticity relations (a a a a ') = (cr_ a a_ a /). There is a simple 
proof of the first relation in ( [2.36| ) based on Fredholm determinants. For the particular 
correlator ( 01/201/2)5 from ( |2.22| )( p723| ) one finds that = K |_ = V, so that 



Det + K (1/2 ' 1/2) ) = Det ^ \) = Det(l + V) Det(l - V) = t+i 



(2.37) 



The second relation in (|2.36| ) is more non-trivial, since it relates a Fredholm determinant 
to a sum of such determinants. We have not found a simple proof of it, however we have 
verified it to lowest non-trivial order. 

In the next section, we will derive differential equations satisfied by the logarithm of 
the SGff correlation functions. For (01/201/2)5 these differential equations can obviously 
be derived from the known Ising differential equation, since log(r-|_r_) = logr_|_ + logr_. 
However, for (01/20-1/2)5 it is not evident how to obtain our result from the known Ising 
results. 

The methods we develop in the next sections can be used to provide a new proof of the 
known differential equations which govern the Ising correlation functions; this is described 
in appendix D. 



3. Differential Equations for Correlation Functions. 

In this section we prove the result stated in the introduction by using the Fredholm 
determinant expressions for the correlation functions derived in the last section. The 
method we follow is a generalization of methods developed in @@]. As will be clear in 
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the sequel, the proof does not rely on the detailed properties of the kernel K, such as its 
a dependence, but rather only on its most salient features. 

The z, ~z dependence of all the quantities we will introduce are a consequence of the 
z, z dependence of the functions e(u) and e(u). The latter are very simple: 



o z e(u) = — 2~e(u), d z e{u) = - — e[u) 
dze(u) = -—e'(u), d^e{u) = -ye'(«), 

where for convenience we have defined: 



(3.1) 



e'(«)^, e'(u) S M. (3.2) 
u u 



Expressions for derivatives of K follow from eqs. ( |3.1|) : 

Tft Tfl 

d z K_ + (u,v) = -—e{u) e(v), d z K + _(u,v) = ~—e(u) e(v) 

Tfl Tfl V / 

d^K_ + (u,v) = -—e'(u) e'(v), d I K + _{u 1 v) = - — e'(u) e'(v). 

Thus, one sees that <9 Z K, djK. are projectors. 
3.1 The resolvent 

We will need the resolvent R for the kernel K, defined by the matrix equation 

(1-R)(1+K) = 1. (3.4) 

As for K, the resolvent R is a 2 x 2 matrix of integral operators. 

For simplicity, consider first scalar kernels of the following form: 

W".»)= £ ° e<,( " )e>) (3-5) 
u±v 

The index (±) refers to the nature of the simple pole of Vt±y Kernels of this type form a 
group with a graded multiplication law. Let [Vr±\] denote the class of kernels of the 
type Vr±). Then using the identity 



11 1/1 1 



u + w w + v u—v\w+v u+w 
one has the following multiplication law: 



(3.6) 



[V ( _ } ] x [V ( _ } ] = [V ( _ } ], [V (+) ] x [V(_)] = [%)], [%)] x [V (+) ] = [V ( _)]. (3.7) 
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Scalar kernels are multiplied in the conventional manner: 

'v<UvW){u,v)= r—VW(u,w)V&(w,v). (3.1 

/ Vn W 



For example 



TO ) («,«) = * • < 3 - 9 ) 



u — V 

This generalizes the study of ||[f|, where only kernels of type V(_) were considered, which 
by themselves form a group. This group property is important since it ensures that the 
resolvent is also an element of the group. In Appendix D, further properties of these kernels 
and their resolvents are presented. 

Let us return now to the construction of R. For future convenience, consider the 
slightly more general problem of finding the resolvent H± for the kernel ±K: 

(1-R±)(1±K) = 1, (3.10) 

where R = R+. Expressing 

**={L If) < 3 - n » 

where D, E, F, G are scalar kernels, equation (|3.10| ) is equivalent to the system of equa- 
tions 

F + DK+- = 0, K-+ =D + FK-+ 

(3.12) 

G + EK- + = 0, K+-=E + GK+-. 

Due to the fact that K = K*, one has F* = F, D l = E, G* = G, where t denotes 

transpose of matrices and kernels. Since and K |_ are of type [V(+)], one deduces 

that F,Ge [V ( _)] and D, E e [V {+) ]. 

The equations ( |3.12|) can be used to express R in terms of K. We first state the result 
then describe how to prove it. Define the functions where 

|+±|-)=2(1-R±) (|), (3.13) 
so that 

£+ = 2(l-F)e, C = -2De 

(3.14) 

= 2(1 - G)e, C = -2Ee. 
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The £ are functions of 2, and w. Functions are multiplied by integral operators in 
the same way as in (|3.8j), e.g. 



dv 



F(u, v)e(v). 



(Fe)(u) 

Then R can be expressed in terms of these functions: 

n( v e+(«)?-(«)-r(«)r(«) 

F(«,v) 
G(u, v) 



(3.15) 









c«)r(«) 


4(u + v] 








4(w — u ) 




f + («)r(«)-r 





(3.16) 



4(u - i?) 

For future reference we will also need 

2e = t + CK+-, = r + t + K+- 
2e = l + + CK-+, = f- + 

These relations are the component form of the matrix relation 



(3.17) 



(3.18) 



the latter being a simple consequence of (|3.13|) . 

The proof that the resolvent takes the form indicated in Q3.16 ) is by direct computa- 
tion. We illustrate this by proving the first relation in ( 3.12 ), F + DK^ = 0, the other 

relations being proved similarly. Using the expression for D in (|3.16|) , one obtains 



(DK+_)(u,v) 



Using now the identity 



1 f dw 



1 



w \u + w 

this becomes 



1 



W + V 



e(w)e(v). 

(3.19) 



u — v 



-(De)(u)e(v) + ^+(u)g+K + _)(v) - \c («)(fJf + _)(«) 



(3.20) 

The relations (|3.14| ) and ( |3.17| ) can now be used to show that the RHS of ( |3.20j ) equals 
—F. 
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3.2 The potentials 

We now have to define certain functions which are essential for the proof, since they 
are directly related to the derivatives of the logarithm of the correlation functions. Define 
the "potentials" 

S± = (e,^ ± )-(lTl), B ± = (e,i ± )-(l T l) 

(3.21) 

C± = (e', e ± > - (1 =F 1), C± = (e - (1 =F 1) 
where ( , ) denotes the scalar product defined for any functions a(u) and &(w) by 

(a,b}= [ —a(u)b(u). (3.22) 
Jo M 

The potentials are functions only of z, z. They are not all independent due to the symmetry 
properties of K. We have the following algebraic relations : 



C-C- - c+d+ = 4 



(3.23) 



The first relation is easily proved using D l = E: 

B_ = -2(e, De) = -2(e, De)* = -2(e, Ee) = S_ . (3.24) 

The second relation is proven in appendix B. 

Let E(#,z) denote the logarithm of the Fredholm determinant : 

S(z, z) = log Det(l + K) (3.25) 

Using log Det(l + K) = Tr log(l + K) , its derivative with respect to z is expressed in terms 
of the resolvent : 

d z E(z, z) = Tt [(1 - K)d z K] (3.26) 



The trace can easily be computed since (<9 Z K) is a projector. From (|3.3|) , (|3.14|) one 
obtains: 

a z E(z, z) = ~ (B_ + B_ + 4) = ~{B. + 2) (3.27) 

To find a differential equation obeyed by E, one next has to derive differential equations 
for the potentials. This will be done in the next sections. 
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3.3 The auxiliary linear system 

Non-linear differential equations for the potentials can be derived by introducing an 
auxiliary linear system of differential equations. This system is obtained by computing the 
derivatives of the two dimensional vectors £ and £ with components ^ and ^ : 

««>=(£(«>)■ MR:!)- (3 - 28) 

The z, z derivatives of these vector functions are computed using their explicit definition 
in terms of the components D, E, F, G of the resolvent ( |3.16|) and the formula (|3.1|) for 
the derivatives of the functions e(u) and e(u). One finds that : 



with 




(3.29) 



2 \B + -u 

mfC-d-+C+d + -2C-d + 



(3.30) 



8u\ 2C+C- -C-C--C+C+ 



The matrices A z and Ay have the same expressions as A z and A- but with "hatted" and 
"unhatted" potentials interchanged. The proof of these equations is given in Appendix C. 

3.4 The differential equations 

From the linear system ( |3.29| ), one can deduce differential equations satisfied by the 



potentials. The result is the following: 



d z C+ = 




d z C. = 




d z C+ = 




d z d- = 




&zB+ = 




&zB+ = 




c\B.= 






-4) = "i C + C + 



(3.31) 



Since all of the equations in (|3.31|) are proven similarly, we illustrate the proof by 



example. From the definition of B + , one has 

d,B + = (d z e^+) + (e,d^+). (3.32) 
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Substituting (|3.1| ) and <9j£ + from ( |3.29|) , one finds 



d^B + = -|<e',£+) - j(C+C + + C_CL)<e, ^) + jC + C.(e, (3.33) 



Noting that (e, £ /if) = (e ; , ^ ± ), and using the definitions ( |3.21|) , one obtains the result in 
( |3.31|) . Consider next d z C + : 

d z C+ = (d z e\t) + (e\d z t) 

m + m ,, . m . (3.34) 

Noting (e',ut; + ) = (e,£ + ), one obtains the desired result. 

One can also find differential equations for the potentials from the zero-curvature 
condition for the connection A z , A^ : 

[d z + A z , dj + A z ] = 0, (3.35) 

and similary for the connection A z , A z . These zero curvature conditions yield a set of 
differential equations which is in fact weaker than, but compatible with, the set of equations 
( p.31| ) . Since we will not need the equations which follow from the zero-curvature condition, 
we do not present the details. 

The last equation in ( |3.31| ), along with (|3.27|) , can be used to express the derivatives 
of the logarithm of the correlation functions only in terms of the C potentials: 

2 

777 . — 

ckd z Z(z,z) = -—C+C+. (3.36) 

Notice that the B potentials can be eliminated from eqs. ( |3.31|) . This leads to differ- 
ential equations involving only the C potentials : 



C-d z C- = C+d z C+ , C+d z C+ = C-d z C- 

9 

777 ^ 

C+(d z ckC.) - (d z C.)(chC + ) = — C-C + (C + C+) (3.37) 

9 

777 ^ 

C-(d z d I C+) - (d z C + ){ckC.) = — C_C+(C_C_) 

Two additional equations corresponding to interchanging C± with C± are derived similarly. 
These equations are similar but not identical to the Hirota equation for the s/2 affine 
hierarchy. 
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In order to solve ( |3.31| ) and (|3.37|) , we introduce a parameterization of the C±,C± 
potentials which solves the algebraic constraint (|3.23| ) : 

C_ = 2e a cosh ip, C + = 2e b sinhv? 

(3.38) 

C_ = 2e~ a cosh^, C+ = 2e~ b sinhp 

with a, b and p auxiliary functions, depending on z, ~z. Inserting this parameterization 
( p.38|) into the differential equations gives the following. The first two equations in ( |3.37Q 
give 

d z a = — tanh 2 p d z b. (3.39) 
Using this equation and its &j derivative the second 2 equations can be simplified to 

(d z dz-a) coth^ — (d z dzb) tanh p — 2d z p> dz(b — a) = 

m 2 (3.40) 
d z djp> = — sinh2<£> — tanhy? d z b d^(b — a). 

The function b can be deduced using Lorentz invariance. Let z = re l9 /2, z = re~ lB /2, 
and consider shifts of 6 by 7. The functions e, e satisfy 



e(e i7 z, e~ ll z, u) = e^ 1 ^' ~ a) / 2 e(z, z, e^u) 
e(e^z, e'^z, u) = e'^+^^^eiz, z, e ll u). 



(3.41) 



From the definition ( |3.21| ) of C+, C+, and making the change of variables u — > e 11 u, one 
finds 

C + = e 2l(Q - Q ' )e /(r), C+ = e" 2l(a - a ' )e /(r), (3.42) 
for some scalar functions /, /. Then, using 

e(z,z,u)=ue(z,z,l/u), (3.43) 

one can show / = / by making the change of variables u — > 1/u. Thus, 

e 2b = £+ = e 4 l( a-a> (3.44) 

C+ 

and 

b = (a - a') log (|) . (3.45) 
Inserting this b into ( |3.39|) and taking the complex conjugate, one deduces 

dz-a = tanh 2 <p d^b. (3.46) 
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The function (p is only a function of r. Thus (|3.40 ) can be written as 



1 „ \ m z . . _ 4(a — a') 



l\2 



d r H — d r ip = — sinh 2(p + 
r I 2 



tanho?(l — tanh 2 ip). 



(3.47) 



Finally, using the equation (|3.36|) , and also (|3.38|) , one obtains 



1 \ TTl 

<9 2 H — d r \ E(r) = — m 2 sinh ip = (1 — cosh 2ip) . 



(3.48) 



This is the result announced in the introduction. Notice that <9 2 <9j£ is only parameterized 
by a single function <p(r), and the differential equation for <p involves only ip itself. 



4. Concluding Remarks 

Though the proof is somewhat lengthy, the main result stated in the introduction is 
strikingly simple. What about away from the free fermion point? Certainly, our derivation 
relies essentially on the fact that the particles in the form factor sum are free. This fact 
is ultimately responsible for the properties of the form factors that allowed us to express 
the form factor sum as a Fredholm determinant. Nevertheless, the simplicity of the final 
result suggests it is worthwhile to investigate whether other derivations of it are possible 
that make no reference to the free fermions. For instance, the SGff theory possesses an 
exact si (2) affine symmetry ||1 5|| , which is known to become q deformed away from the free 
fermion point It would be very interesting to be able to derive the result presented 
in this paper solely from the si (2) symmetry, as such a derivation may be amenable to 
(/-deformation. 
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Appendix A. Form factors from monodromy. 

Here we describe how the form factors in section 2 can be derived from braiding 
relations with the fermions. As in the ising case||, the braiding relations induce a Bo- 
goliubov transformation on the creation-annihilation operators in momentum space. The 
field, which is the unitary operator implementing this transformation, is therefore the 
exponential of a bilinear form in the fermions. 

Let = ( ^ ± ^ be the Dirac field satisfying the Dirac equation: 

d z ifj± = irmfj±, dgij)± = —imtyj^ (A.l) 
Its mode decomposition in momentum space is: 

V>±(aM) = ±V^ J d6e e/2 (c ± {9)e- ip - x - cj(6»)e ip - x ) 

$ ± (x,t) = TiV^ J dOe' / 2 (c ± {d) e - ip - x + c£(0)e ip ' 

with p ■ x = m[x sinh 6 — t cosh 9) . 

Let a a =: exp (zct$(0, 0)) :. This field is non-local with respect to the fermions, and 
therefore has non-trivial equal time braiding relations with them: 

a a *±(y) a" 1 = *±(y), for y < 

(A.3) 

a a ¥±(y) a- 1 = e ±l2 ™ *±(y), for y > 

Here, the fermions &±(y) are at t = 0. 

Equation ( |A.3|) induces a Bogoliubov transformation on the C±(9Ys and C^(6*)'s. To 
describe this transformation, one has to Fourier transform a a Sfr±(y)a~ ■ Decomposing it 
into the sum of two integrals, one from — oo to zero, the other from zero to +oo, and using 
the braiding relation (|A.3|) , we get: 

= ! %* k y* ± (y) + (e ±i2 ™ - 1) / p-e*V9 ± (v) 
J-oo 2tt Jo 27r 
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Paramatrizing k as k = m sinh f3, the integrals are easily done. We find: 



{a a C ± ((3) a- 1 ) = C ± ((3) + 



,±i27TCK 



Am 



+00 



d6 



C±{9) 



sinh( 



2 



4(g) 

cosh(^ 



(A.5) 



t 



Here e — > + . A similar relation holds for C_J_(/3) since cr Q is a unitary operator. 

Equation (|A.5|) is a linear transformation of the creation-annihilation operators which 
preserves their commutation relations, since it is defined by conjugation. Therefore it is 
a Bogoliubov transformation. We need a few basic facts about these transformations, cf. 
[2TJ. Let Ci and c\ be fermionic operators with {Ci,cj} = Sij. Consider the linear 
transformation C{ — > specified by: 

j 

This is a canonical transformation if the operators Di and d\ also satisfy the fermionic 
relations {-Di,-Dj} = 8^. This amounts to a few constraints on the matrices X and Y. 
The unitary operator U which implements this transformation, i.e. Di = UCiU -1 , is given 
by: 



(A.6) 



where Af 



Det(xtx) 



If : ex P \l E {Cj^id + zdjRjiG + CjF^) j : 



(A.7) 



-1/4 



The unitary operator ( |A.7|) exists whenever this normali- 

(A.8) 



sation constant is finite. In eq. ( |A.7 ), the matrices F c , R and F are defined by: 
F C = Y(1 + R), X(l + R) = l, F=(1 + R)Y 



These matrices are the two-particle form factors of the operator U, e.g. 



1 



(A.9) 



Therefore, the operator U is completely characterized by its two-particle form factors. 

In our case, the fermion operators C v {6) have a discrete index v = ± and a continuous 
index 9 taking values from —00 to +00. The matrices X and Y are therefore kernels. From 
eq. ( |A.5| ), we deduce that: 

X^' (9, 9') = 8 V >"' (5(9 - 9') + V u (d, 9')) 

Y u,v '(9, 9') =6^- v 'Y v (9,9') 



(A.10) 
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with 

/ p iv2-Ka 1 \ I 

v u {ex ' 



Y u (9,9') 



4%i J smh( e - e '- i£ ) 

x (A.ll) 

iv lira 1 \ 1 



4tu y cosh(^) 

To find the expression of the "spin" field a a , we first have to find the resolvent R. It 
is of the form 

H v,u ' (6, 9') = R v {6, 9') (A.12) 
The kernel R v is the resolvent of V v : 

(1 + V V )(1 + R V ) = 1 (A.13) 

We claim that a solution to ( |A.13|) is : 

/ 1 _ p iv2-Kot\ -1,(0-0') 

^'•^ = (-6S-) (A ' 14) 

To prove this, we have to compute the product (V U R U ): 

/oo 
d/3V v (9,P)R v (P,9'). (A.15) 
-oo 

This involves an integral over the intermediate rapidity (3 from — oo to +oo. Since the 
function to integrate is (z27r)-periodic in f3 up to a phase, which is equal to e _l!y27ra , we 
can compute its integral by considering the contour in the f3 complex plane to be the sum 
of 2 contours, one along the real axis at Im(/3) = 0, the other along the line Im(/3) = — 2ni. 
The integral is then given by the residue theorem; the poles, inside the contour, which are 
all simple poles, are at j3 = 9 — ie and at j3 = 9' — is. Suming the residue contributions 
gives: 

V u Rv + V u + R u = (A.16) 

as claimed. 

The matrix kernel F has the following form: 

F v ' v '(6,6') = 8~ u ' v ' F u (6,6') (A.17) 

with F v = (1 + Rv)Y v . They are computed in a similar way using contour integrals. We 
find 

F .M = V <E^<-?£L (A.18) 
2,171 cosh^^-^-j 

This coincides with the two particle form factors given in the text. The crossed form 
factors F c are derived in the same way. 
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Appendix B. Relations between the potentials. 

Let the functions e',e' of z,z and u be as defined in (|3.1|) . Analagously to ( |3.14|) , 
define additional vector functions of z,~z and u as follows: 

(+ = 2(1 - F)e', C = ~2De' 

(B.l) 

C+ = 2(1-G)e', C = ~2Ee'. 

These functions can be linearly related to the ^ functions, the coefficients being the 
potentials C±,C±. The result is 

r(w) = ^(c_r(«)-c f e + («)) 

( + (u) = ±{-C^+(u) + C + C(u)) 

2u (B.2) 

r(«) = -(c_r( w )-c + ^ 

2u V 



c + («) = ^pi + w + c + r(«)). 

All of the above relations are proven similarly, so we illustrate the technique with one 
of them. From the definition ( |B.lj ), and the explicit form of D in ( |3.16| ), one has 

C(«) =-\f— (-^-) (* + («)£ + («0 ~ C(u)C(w)) e(w). (B.3) 
2 J w \u + w w J V / 



Using the relation ( |3.6|) with v = 0, 

C" («) = ^ (4(i>e)(«) + r + 2) - £+(«)£+) . (B.4) 

Finally, using the expression for De from ( p.!4|) , one obtains the first equation in (|B . 2| ) . 

The relations (|B.2|) can be used to prove the constraint (|3.23|) amoung the C- 
potentials. Namely, using the definition ( p.l| ) for (~ , and the expression (|B.2|) for (~ 
in terms of £ , one has 

2(e, J De') = -(e,^-(C_r-C i + e + )) 
lu 



-i(e_<e',r>-^+<e , ,e + >) (B.5) 
1 -(p-{C-+2)-d + C + ). 



On the other hand, 

2(e, De') = 2(e, De')* = 2(e ', £e) 
= -(£_+ 2). 

Comparing (|B.5| ) and ( |B.6| ), one obtains the constraint ( |3.23| ). 
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:b. 6 ) 



Appendix C. Proof of the linear system. 

In this appendix, we prove the equations (|3.29|) , (|3.30|) . Multiplying (|3.13|) by 1 ± K, 
one obtains 



(1±K)( )■ (CI) 



Taking the ^-derivative of this equation, then multiplying on the left by 1 — R±, leads to 
9 > ( ) = 2d - «W (tl) =F (1 - R ± )S,K ( £ ± C j . (c.2) 

From (|3.3|) and ( |3.14j ), one has 

((i-R ± )a)(^)-^^ +(#) ±rwe(,)J- (Q3) 

Performing a simple multiplication gives 

(1 - R ± )d z K ( g± C ) = ± 5_ ± 2) ) - ± B_ ± 2) ( C_ ) . 

(C.4) 

The other term on the RHS of ( |C.2| ) may be simplified to 

ri jyJd z e\_ m f (l-F)(ue)TD(ue) \ 

( " ±J Va z e J " " T V T-B(we) + (1 - G)(«e) J ' { ^ b) 

Expressions of this kind can all be computed similarly. This requires commuting the 
multiplication by u and the operators E, D : F, G. Consider for example 



/div 
— Div, w)(w + v — v)e(w) 
w 



-v(De)(v) + I — (w + v )D(v , w)e(w) 
w 



(C.6) 



Inserting the explicit expression (|3.16| ) for D in the second term, we obtain 

(D(ue)) (v) = V -C(v) + \ (t(v)B + - C(v)(B- + 2)) • (C.7) 



Similarly, 



(E(ue)) (v) = V -C(v) + \ (t + (v)B + - C(v)(B_ + 2)) 
((1 - F){ue)) (v) = ~t+(v) + i {£ + (v)(B_ + 2) - £~(v)B + ) (C.8) 
((1 - G)(ue)) (v) = ~£+(v) + \ (j + (v)(B_ + 2) - C(v)B + ) . 
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Collecting all terms on the RHS of ( |U.2j ), and solving for d z ^, yields the d z part of 

the result ( ^29|) . 

Consider now the z derivatives. By the same reasoning as above, we begin with ( |C.2| ) 
but with d z replaced with <9j. One needs 

where we have used the definition (|B . 1|) . As in ( |C.4| ), one also has 
(1 - R ± )^K f £±C ) = -^(C + ± C. ± 2) (±f ) - ^(C + ±d.±2) (£_ ) . 



Putting all of this together, one finds that the result simplifies greatly: 



(CIO) 



(C.ll) 



Finally, using the expressions (|B.2|) for C ± 5C ± i n terms of ^ ± ,^ ± , one obtains the rest of 

0, dpo|) . 



Appendix D. The Ising case. 

In this appendix, we provide a new proof of the differential equation for the Ising 
correlation functions ( |2.32| ), and give more details concerning the properties of the kernels 
of the type ( |3.5| ). These correlators are expressed as the Fredholm determinant for the 
scalar kernel ( |2.35| ) which is in the class of integral operators (|3.5| ). Rather than repeat 
what was done above for the sine-Gordon theory, we present the proof in the case where 
the kernel is replaced by a finite dimensional matrix. More specifically, we first cutoff the 
integrals by replacing J °° du with J Q A du for A finite. Consider then discretizing the u 
variables into N intervals, where in each interval u takes the value Wj, i = 1, .., N : 

/ dufiu) - A«X)/(«i), (D.l) 

^0 i=l 

where u n = nAu, and NAu = A. The integral operators arise in the continuum limit 
where — > oo, Aw — > 0, with NAu kept fixed; indeed, this is the way in which the 
Fredholm theory is usually developed. In the end, we take A to infinity. We do this in 
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order to illustrate a basic point: the above proofs do not rely on being in the continuum 
limit but are actually valid for finite N. Of course, all of the proofs below can be done 
also in the continuous case as was done above for sine-Gordon. 
Therefore, we consider two tau functions, 

r± = det (1 ± V) , (D.2) 

where V is an iV x N matrix with entries, 

Vij = (D.3) 

Ui + Uj 

with (ui) a set of complex numbers. The (e^) are functions depending on the parameters 
Ui as well as on the coordinates (z,z). As explained above, the actual continuum Fred- 



holm determinant ( 2.33 ) is obtained by letting N go to infinity, identifying the u^ as a 
discretization of the continuum variable u, and replacing the infinite sums with integrals. 
We assume that the (z,~z) dependence of V comes only from the functions through the 
following formula : 

_ muj 

_ A_ (D.4) 
2m 

It is well known that this implies that r± are the tau functions of the N soliton solutions 
of the classical sine-Gordon equations. This is the result that we will rederive using our 
technique. The method is similar to the case described in the main text but all the 
computations are actually simpler since the kernel is a scalar. The key point is the group 
property ( |3.7|) with the 21 2 graded multiplication law. 
We first need the resolvent R± of the matrix (±V) : 

(1-R ± )(1±V) = 1. (D.5) 

Similarly to the integral operators ( |3.5| ), the matrices of the form 



form a group with the graded multiplication law sketched in eq. (|3.7|) . Therefore, the 
resolvent is an element of the group. From the equation (p.5|) , it is clear that R± has 
components of both type [V(±)]. Namely, we have 

R± = H±F (D.7) 
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where H is matrix of type [Vj_)] and F a matrix of type . Since the matrix ( p.3| ) is 

of type [V( + )], eq. (|L).5|) is equivalent to : 



VF + H = 0, VH + F = V. (D.8) 

These equations can be solved, giving the following expression for the matrices F and H : 

f-f + - f+f- 

%J 2(m — uA 

+ _ _ + (D.9) 

13 2(ui + Uj) 

with 

ft = [(1 - R±)el . (D.10) 



The proof consists in checking directly the relations (|D.8| ) using ( |D.10| ). 



Besides the vector (e^) and (ft), we also define vectors (e[) and (f i ) by : 

f'f = [(l-R±)el, 4=X (D.ll) 
We need to introduce the "potentials" . They are defined by : 

P± = (eJ ± ), Q± = (e'j ± ), (D.12) 
where the scalar product of two vectors with components (oj) and is : 

(a,b) = J2 a i b i- (D.13) 

i 

The /' ^ and ft are not independent, but are related by the potentials. We have : 

f't = -(lTQ±) ft- (D.14) 

Ui 

The proof goes as follows. By the definition ( p.ll| ) we have : 

f't = ^-[(H±F)e / ] i . (D.15) 
Then, using the formula (p.9| ) for H and F, we deduce that : 



(He'), = ±- (f~Q + - /+Q_) + 1 (#e), 
(F e ') i = ^ : (rQ + +/ + Q-)-^(F e ) 4 . 
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(D.16) 



Substituting the latter into ( p.!5|) proves (|D.14| ). 

The relation ( p.!4j) implies that the potentials Q± are not independent. Indeed, 
taking the scalar product of /' ± with (e^) and using (e, /' ± ) = (e', f^} = Q±, we obtain : 



;i-q+)(i + q_) = i. 



(D.17) 



The derivatives of the logarithm of the tau functions ( p.2| ) are given by the potentials. 
Using d z logr± = ±tr [(1 — R±)d z V] and ( P-4j ), we obtain : 



d z logr± = T y(e,/ ± ) 
^logri = Ty(e',f ) 



(D.18) 



The next step consists in finding the auxiliary linear system. This is a system of linear 

f f+\ 

differential equations for the N two-dimensional vectors [j- Let us first describe the 
result and then give the proof. The system is : 

' f+\ m f -P+ - P~ 



d z 



fi 



Ua 



+ 



(D.19) 



1 fr 



rr^f (l_g+)2W/+ 

2 U A( 1 +Q") 2 )\fr 



Let us first prove the d z equations. Taking the derivative of ( p.l(J| ) and using the expression 
( p.9| ) for the resolvent, we obtain : 



in 



[(1-R±)(ue) 



(D.20) 



Now, the last term in ( p.20| ) can be simplified using ( p.9| ) 



[(1 - RtXveftt = TftP^ + «i [(1 " R T >] 



(D.21) 



Inserting this relation into (|D.20|) gives the d z part of (|D.19| ). The <9j equations follow 
from : 

&zft = -y (1 =F Q±) ft, (D.22) 



which is a direct consequence of (p.4|) and of the definition ( p.ld|) . Expressing (ft) via 
equation ( p.!4| ) gives the dz part of (|D.19| ) . 
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We now relate the potentials P± to the Q±. The d z derivative of Q± is 



d z Q± = (d z e , J ± } + (e',d z f ± ). (D.23) 



Using (|D.4|) and the linear system ( p.!9|) , we obtain : 



777 

d z Q± = -— (1 =F Q±) (P+ + P-). (D.24) 
Solving the algebraic relation ( |L) . 1 7| ) by parametrizing the potentials Q± by : 

1TQ± =exp(T^), (D.25) 
for some auxiliary function (p, we deduce that : 

ID 

-( P+ + P_) = -d z <p. (D.26) 



Inserting this parametrization in (p. 19 ), one recognizes the Lax pair of the sinh- 



Gordon model. Therefore, its zero curvature condition gives : 

9 

Tfl 

dxdjgip = — sinh2</?. (D.27) 

The last step consists now in relating the tau functions ( p.2|) to the function ip. Using 
([DT8D , one finds that : 

TTl 

9 z ajlog(r + r_) = — (1 - cosh 2<p) 

fr + \ m 2 (D - 28) 
d z dj\og ( — J = — sinh2v?. 



These are the equations of motion of the affine s/2 Toda theory p2 |. 



Note that since the Ising correlators are r± themselves rather than their product 
and quotient, the above parametrization of the Ising correlators is somewhat awkward 
in comparison with the sine-Gordon case. This suggests that the differential equations 
perhaps have a more natural explanation in the sine-Gordon case. 
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